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Abstract 

We have compared the intensity distribution of molec- 
ular line emission with that of dust continuum emission, 
and modeled molecular line profiles in three different pre- 
protostcllar cores in order to test how dynamical evolu- 
tion is related to chemical evolution, and whether we can 
use different chemical tracers to identify specific dynami- 
cal evolutionary stages. We used dust continuum emission 
to obtain the input density and temperature structures 
by calculating radiative transfer of dust emission. Our 
results show that chemical evolution is dependent on dy- 
namical processes, which can give different evolutionary 
timescales, as well as the density structure of the core. 
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1. Introduction and Observations 

Low mass star formation has received more detailed study, 
both in theory and observation, compared to massive star 
formation. An evolutionary seqence of low mass star for- 
mation based on the spectral energy distribution (SED) 
has been developed (Adams et al. 1987). However, the 
initial conditions of collapse are still poorly understood. 
Preprotostellar cores (PPCs) are believed to be gravita- 
tional bound, but they have no central hydrostatic proto- 
star (Ward-Thompson et al. 1994). Therefore, PPCs can 
be considered as the potential sites of future star forma- 
tion and give us chances to probe the initial conditions 
of star formation. Evans et al. (2001) modeled dust con- 
tinuum in three PPCs using Bonnor-Ebert spheres and 
showed that the dust temperature decreases toward the 
center because PPCs are heated only by the interstellar 
radiation field, and the radiation from the cores is opti- 
cally thin. Therefore, they argued that the actual density 
structure of a preprotostellar core is more centrally con- 
densed compared to the intensity distribution of the dust 
continuum emission. 



Several studies on chemical evolution in PPCs have 
been done by Bergin & Langer (1997), Aikawa et al. 
(2001), Caselli et al. (2002), and Li et al. (2002). They 
tested different dynamical scenarios, but all included the 
interaction of gas phase molecules with the surfaces of 
dust grains as well as gas phase chemistry and showed 
that some molecules such as CS and CO are depleted 
substantially from the gas phase as the density increases. 
Those molecules have high binding energy onto the sur- 
faces of dust grains, so their timescales for evaporation 
by thermal energy or cosmic rays are much longer than 
those of molecules that have lower binding energy. For 
instance, nitrogen bearing molecules such as N 2 H + and 
NH3, which become abundant in later stages of dynam- 
ical evolution and have low binding energy, do not show 
as significant depletion as CO and CS. Comparison of the 
abundance change of a molecule with time and the dif- 
ferential distributions of molecules in a given time seems 
to be able to trace the dynamical evolutionary stages. 
We tested this idea by comparing the line emission from 
several molecules that have different chemical evolution 
timescales with dust emission, which can trace physical 
conditions, in three PPCs: L1512, L1544, and L1689B. 

We observed the three PPCs in C 18 (J = 2 - 1 and 
J = 3 - 2), C 17 (J = 2- 1), DCO+ (J = 3 - 2), HCO+ 
(J = 3-2), and H 13 CO+ (J = 3 — 2) with the 10.4 m tele- 
scope of the Caltech Submillimeter Observatory (CSO) at 
Mauna Kea, Hawaii from 1995 to 2002. We also observed 
these cores in H 13 CO+ (J = 1 - 0), N 2 H+ (J = 1 - 0), 
and CCS (Nj = 43 — 3a) with the 45 m telescope of the 
Nobeyama Radio Observatory in Japan on January 2002. 
See Lee et al. (2003) for details about observations. 

2. Simple Analysis 

We compare, first, the H 2 column densities calculated 
from molecular lines (C 18 J = 2- 1 and C 17 J = 2-1) 
and dust continuum, which can be used as the tracer of 
column density least affected by chemistry (Figure l.a). 
We assume that temperature and abundance are constant 
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along the line of sight, that gas and dust have the same 
temperature, that all levels are in LTE, and that lines are 
optically thin. The H2 column densities calculated from 
C 18 J = 2 — 1 are smaller than those calculated from 
850 /jm dust emission by a factor of 10 to 25. Also, the 
C 18 emission shows no peak at the position of the peak 
of dust continuum emission. If this line is optically thin, 
so that it traces all material along the line of sight, the 
difference would be caused by the depiction of CO. The 
optical depth of C 17 J = 2 — 1 can be calculated by 
fitting the relative strength of its hyperfine components, 
and the calculated optical depths in three PPCs show 
that the optical depth ofC 17 J = 2 — 1 is negligible in 
all three cores, but C 18 J = 2 — 1 is optically thick in 
L1544 and L1689B (see §4.2.1 of Lee et al. 2003). There- 
fore, if we use an optically thick line, the depletion factor 
can be overestimated. The comparison between the H2 
column densities calculated from C 17 J = 2 — 1 and 850 
/im dust emission shows that the CO molecule is signifi- 
cantly depleted in L1512 and L1544 but not in L1689B. 

We also compare the distributions of the intensity of 
molecular emission and dust continuum emission (Figure 
Lb). L1512 is the weakest, and L1544 and L1689B are 
similar in dust continuum. However, the CCS intensity 
of L1689B is more similar to that of L1512. Even though 
the distribution of CCS is very different from that of the 
dust continuum, it does not show the central dip expected 
if the molecule is depleted significantly as some chemical 
modelings show. We suggest that this trend would be 
caused by the significant depletion of CO, helping C + to 
react more with molecules that are not oxygen-bearing 
species (Li ct al. 2002). The intensity of N 2 H+ follows 
dust emission well in L1544 and L1689B, but it is weaker 
in L1689B than in L1544. The difference is possible be- 
cause CO, which is a destroyer of N2H + , is much more 
depleted in L1544. On the other hand, in the central re- 
gion of L1512, the N 2 H+ distribution is flatter than the 
dust emission indicating the possible depletion of N 2 H+. 

3. Detailed Models 

The simple method used in the previous section has limi- 
tations for quantitative analysis because it cannot include 
the variation of the excitation temperature and the abun- 
dance of a molecule along the line of sight. Therefore, we 
used the Monte Carlo (MC) code and the virtual tele- 
scope simulation code developed by Choi et al. (1995) 
to calculate radiative transfer of molecular lines and to 
simulate specific molecular line profiles, respectively. For 
the simulation, we adopted two physical models: Bonnor- 
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Figure 1. (a) compares the H2 column densities calculated 
from the 850 [im dust continuum emission, C 18 J = 2 — 
1 (without the correction ofr), and C 17 J = 2 — 1 (with 
correction of t) through the cuts marked in Figure 2 of 
Lee et al. (2003). (b) compares the integrated intensities 
of CCS Nj = 4 3 - 3 2 and N 2 H+ J = 1 - lines with S 850 
through the same cuts in (a). Here, Ss5o is shifted by 0.9. 



Ebert spheres (Bonnor 1956, Ebert 1955), which give only 
density structure; and Plummer-like spheres (Whitworth 
& Ward-Thompson 2001), which give density and veloc- 
ity structures and can be used for optically thick lines 
to simulate their self-absorbed line profiles (see §5.2 of 
Lee et al. (2003) for details). We calculated the radia- 
tive transfer of dust emission to obtain the temperature 
structure (Evans ct al. 2001). For the structure of molec- 
ular abundances, we tested several functional forms; the 
step function, which defines a depletion factor (fo) and 
a depletion radius (r£>), was the best to fit line profiles 
from the center to the outer regions in the three PPCs 
(Fig. 5 in Lee et al. 2003). We modeled C 18 0, H 13 CO+, 
DCO + , and HCO + lines using a step function of each 
abundance in order to test how much those molecules are 
depleted within what radii. Figure 2 shows that the de- 
pletion radius is well constrained around 0.075 pc, but 
the fractional depletion is constrained only to be greater 
than about 25 in the C 18 lines of L1512. This trend oc- 
curs in every molecule that is significantly depleted, so we 
consider fo to be the lower limit of the depletion factor. 
In other words, the data are consistent with complete de- 
pletion inside some Td- Figure 3 shows the best-fit model 
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Figure 2. The distribution of the reduced \ 2 of models that 
have different rp and fo in C ls O 2—1 and 3—2 toward 
L1512. The contour levels are from 5 to 95 with the in- 
terval of 10. 

Table 1. The results of MC models of molecular lines 
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a We could not calculate ro and fo of DCO + in L1512 because 
we have only one spectrum toward the center. 

of C 18 J = 2 - 1 in L1512. The results of modeling in 
the three cores are summarized in Table 1. 

According to our results, CO is substantially depleted 
in L1512 and L1544 but not in L1689B. In addition, the 
depletion radius of CO in L1512 is greater than that 
in L1544. The C 18 lines in L1689B might be affected 
by surrounding warm gas (see §5.3.2 of Lee et al. 2003). 
HCO + is also depleted significantly in L1512 and L1544, 
but the depletion radii are smaller than those of CO. On 
the other hand, DCO + does not show significant deple- 
tion in any of these three cores. We might not detect the 
depletion of DCO + because our observational resolution 
was poor. 

4. Discussion 

Our simple analysis and detailed modeling of molecular 
line profiles combined with dust continuum emission show 
that L1512 is least centrally condensed, but CCS, CO, 
and HCO + are significantly depleted. Even N 2 H+ is pos- 
sibly depleted in L1512. On the other hand, L1689B has 
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FigureS. The best-fit model of C 18 O J = 2 - 1 in L1512. 
The arcseconds marked in panels represent the angular 
distance from the dust peak. 

Table 2. The results of this study 
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a similar density structure to L1544, which is quite cen- 
trally condensed and has substantial depletion of CCS, 
CO, and HCO + , but no molecules except for CCS show 
significant depletion. This result may indicate that L1689B 
has evolved too fast to have time to evolve chemically. 
Therefore, we suggest that the stage of dynamical evolu- 
tion of a core cannot be simply probed by its chemical 
status because chemical evolution depends on the size of 
dust grains (Caselli et al. 2002, in preparation) or the 
absolute dynamical timescale, which is different for each 
dynamical process, as well as its density structure that 
shows the relative dynamical evolutionary stages (Table 
2). 
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